Forced tearing of ductile and brittle thin sheets 
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Tearing a thin sheet by forcing a rigid object through it leads to complex crack morphologies; a single oscil- 
latory crack arises when a tool is driven laterally through a brittle sheet, while two diverging cracks and a series 
of concertina-like folds forms when a tool is forced laterally through a ductile sheet. On the other hand, forcing 
an object perpendicularly through the sheet leads to radial petal-like tears in both ductile and brittle materials. 
To understand these different regimes we use a combination of experiments, simulations and simple theories. 
In particular, we describe the transition from brittle oscillatory tearing via a single crack to ductile concertina 
tearing with two tears by deriving laws that describe the crack paths and wavelength of the concertina folds and 
provide a simple phase diagram for the morphologies in terms of the material properties of the sheet and the 
relative size of the tool. 

PACS numbers: 62.20.mm, 68.60.Bs, 46.15.-x 



The failure of thin sheets by fracturing, tearing and peeling 
occurs naturally on a range of scales - from the everyday world 
of opening an envelope or other package to the grounding of 
ships (3, from the failure of armor plating to the cracking of 
geological structures Q. The study of the fracture and tearing 
of thin sheets is challenging as it couples the geometry of large 
out-of-plane deformations to the failure of the material. De- 
formations leading to in-plane stretching, for example because 
of constraints in the far field, play a crucial role in determining 
the resulting complex crack morphologies. While early work 
focused primarily on the role of plastic deformations (U 
there has been a recent surge of interest EHU in understand- 
ing these questions quantitatively in brittle elastic materials. 
This is driven in large part by the focus on understanding the 
mechanical behavior of thin films at the meso/nano scale, in 
developing of assays for the measurement of material proper- 
ties and creating new functional structures ifTOlfTTl . 

An important class of fracture and failure in thin sheets is 
that due to the forcing of a solid object or tool through it, 
either in the plane of the sheet or transverse to it. Forced 
tearing of a relatively brittle sheet, made of acetate, by the 
motion of a tool in the lateral direction leads to oscillatory 
fracture (Fig. [T^) if the tool diameter D is large in compar- 
ison to the sheet thickness h. For this system, various levels 
of approximation yield a hierarchy of models for crack paths 
lEHZL all of which show that the amplitude and wavelength 
of the oscillation scale with D, consistent with experimental 
observations. Furthermore the dynamics of the single crack 
in this case shows a characteristic stick-slip behavior associ- 
ated with the transitions between bend-dominated and stretch- 
limited motion of the tool. Given the strongly geometric fla- 
vor of all models for this phenomena, a natural question is the 
following: how would a ductile sheet that can deform plasti- 
cally before fracture behave when forced similarly? Here we 
attempt to synthetically understand the phase space of the dif- 
ferent crack morphologies that arise in the tearing of a thin 
sheet as a function of the type of forcing, the geometry of the 
sheet and its material properties. 

In Fig. [TJ) we show the results of fracture in a thin sheet 



of paper driven by the lateral motion of a tool; paper is more 
ductile than acetate, and this leads to disordered version of the 
oscillatory fracture seen in acetate sheets. In Fig. [TJ; we show 
the results of a similar experiment carried out with an alu- 
minum sheet, which is far more ductile, and see two marked 
differences: the sheet fractures by forming two cracks, and 
a series of irreversible periodic folds forms along the crack 
edges. If the same experiment is now carried out in the pres- 
ence of a rigid substrate that supports the thin aluminum sheet 
using a tool that leans forward in the direction of motion, we 
see a beautiful 'concertina' tearing pattern (Fig. |3^) first stud- 
ied in the context of the failure of relatively thick metallic 
plates 1 3 1, and understood using a plastic kinematic deforma- 
tion model. For very thin sheets or for materials with a high 
yield stress, however, this analysis is not valid since the defor- 
mations ahead of the crack tips are predominantly elastic. 

To understand the role of ductility in thin sheet fracture, we 
start with numerical simulations that combine finite elastic- 
plastic deformation with a simple fracture criterion to synthe- 
size the variety of observed fracture morphologies in brittle 
and ductile sheets. Our simulations are based on a discrete- 
element model that consists of mass points arranged in a ran- 
dom lattice with average area h 2 per site, where h is the thick- 
ness of the sheet fl2l . Lattice sites have translational and 
rotational degrees of freedom, and each pair of neighboring 
points is connected by elastic-plastic element with bending, 
shear and tensile stiffness, similar to one used to study crum- 
pling [13 ] and brittle fragmentation fl4l . Lattice sites have 
six neighbors on average. We assume an ideal elastic-plastic 
stress-strain relation with Young's modulus E and yield stress 
(jy . Elements are forced to break when their tensile or shear 
strain exceeds j s . For ductile materials (J y /E < j s , i.e. the 
material deforms plastically before fracture. For cr y /E <C 7 S 
we obtain the toughness T ~ a y h typical for thin metal sheets 
|[T5lL while in brittle case T « Ehj 2 [12]. In our simulations 
T is the sum of energy of broken elements and plastic dissipa- 
tion per unit area of crack extension. The above estimates for 
T provide a convenient way to relate our simulations to exper- 
iments; however the connection of 7 S to real fracture strains 




FIG. 1 : Complex morphologies of cracks in a sheet of (a) brittle poly- 
meric material, (b) paper, and (c) aluminum film torn by a cylindrical 
tool that moves in the plane of the sheet. 



is more complicated. Self- avoidance of the sheet is generated 
via an elastic repulsion with range h between any two non- 
neighboring lattice sites. The dynamics of the system is simu- 
lated by solving Newton's equations of motion. We start with 
a long sheet with dimensions 320h x 1280/i that is clamped 
along its lateral edges with no pre-strain and a small notch is 
introduced on one of the short edges via which a tool is intro- 
duced to cut the sheet. The tool is moved with a given velocity 
that is small enough that the inertia is negligible. We model 
the interaction between the tool and the sheet using a simple 
Coulombic frictional law, with a coefficient fi for both static 
and dynamic friction. 

We first replicate numerically the experiments on cutting 
a sheet with a cylindrical tool of diameter D that ploughs 
through a sheet of thickness h (D > h) laterally l4ll5l[T2ll. 
For a brittle sheet we observe an oscillating crack path as 
shown in Fig. [2^. Following an initial transient the motion 
of the crack settles into that observed in experiments (Fig. 
[T^). Furthermore, the cutting force lH2l is highly oscillatory 
as elastic energy is stored and released in dynamic bursts of 
fracture (U [6) . We confirm that fracture occurs by in-plane 
stretching when the crack tip is at either extreme laterally and 
by out-of-plane shear when the tip is close to the center-line 
13). We use 7 S = 0.4 to obtain Y/Eh ~ 10~\ a value typ- 
ical for polymeric sheets used previously (U|6). The relative 
toughness here is much higher than that of ceramic materials 
and glass, where forced failure occurs through brittle fragmen- 
tation. 

Keeping the geometry of the sheet and tool the same, we 
now vary the onset of plastic yielding by reducing the yield 
strain cr y /E so that it is smaller than the fracture strain 7 S . Os- 
cillatory motion of a single crack persists till cr y /E « 7 s /2; 
as the yield strain is reduced further, the tool propagates by 
forming two parallel cracks rather than a single oscillatory 
one. Our simulations indicate that this transition is indepen- 
dent of D when D ^> h. However, forD^h stable straight 
cracks have been reported for both brittle |4] and ductile (T) 
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FIG. 2: Simulations that reproduce the experiments shown in Fig. 1. 
(a) An oscillating crack in a brittle sheet, (b) Disordered tongues in 
a weakly ductile sheet with a y /E = 7 s /4. (c) Periodic concertina 
tearing in ductile sheet with a y /E = 7 s /40. Contact between the 
tool and sheet has a friction coefficient /x = 0.25 and tool diameter 
D = lOO/i. (d) Phase diagram of tearing morphologies as a func- 
tion of tool size and material behavior characterized by the relative 
magnitude of the yield strain a y /E to the fracture strain 7 S . 



sheets. While the discrete nature of the lattice limits accuracy 
in this limit, we find that for D < lOh crack paths are also 
straight in both the brittle and ductile case. 

For a y /E<j s /2 and D ^> h the nucleation of new cracks 
arises when a single crack is driven away from the middle-line 
by the advancing tool. In a brittle sheet, this is followed by a 
dynamical burst returning the crack tip back to the middle. In 
ductile sheet, however, this step is hindered by plastic defor- 
mation and a new crack nucleates at the point diametrically 
opposite where the crack tip meets the tool. Then, both cracks 
propagate simultaneously while extruding a tongue of mate- 
rial ahead of the tool. Whether or not this tongue folds up 
into a concertina pattern depends on the ductility of the ma- 
terial. For weakly ductile materials, the tongue bends out of 
the plane without forming permanent folds (Fig. [5J)) resulting 
in fracture patterns seen in torn paper (Fig. [lj>). For strongly 
ductile materials, with cr y /E <C 7 S these tongues appear on 
alternating sides of the tool in a form of short concertina folds 
(Fig. [2]:), as seen in torn aluminum films (Fig. [ij). The 
size of the folded regions scales with D and increases with 
increasing friction between the tool and the sheet. Finally, by 
introducing a pair of initial cracks on both sides of the tool 
or a rectangular tool, we find that a stable series of concertina 
folds results. Thus, although the morphology of tearing of 
thin sheets is dependent on tool shape, contact friction and the 
mechanism of crack initiation, at a primitive, coarse level the 
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FIG. 3: (color online). Concertina tearing 
of (a) aluminum film and (b,c) simulated 
ductile sheet on a rigid substrate. In (b) the 
process starts with the formation of a wrin- 
kle of wavelength A that is then stretched 
(c) before the cracks propagate. Color indi- 
cates elastic stretching energy density, (d) 
Experimental and simulated wavelength of 
concertina folds are plotted as a function of 
width w of the tongue and compared with 
X/h ~ (w/h) 2/s . Each plot includes data 
from three (experiments) or two (simula- 
tions) independent measurements. In each 
simulated case a y /E — 7 s /20. (e) Tongue 
width w is plotted as a function of length 
x of the tongue and compared with w/h ~ 
(x/h) 3 / 4 . Simulated x and w are magni- 
fied by factors 8 and 8 3//4 , respectively, for 
better visual comparison. 



patterns are characterized by material behavior; ductile tear- 
ing is characterized by the formation of a tongue between two 
cracks that may fold up ahead of the tool, while brittle tearing 
is characterized by the oscillatory motion of a single crack. 
This simple phase diagram is summarized in Fig. |2|l in terms 
of the tool size and material behavior of the sheet. 

To analyze concertina tearing in more detail, we numeri- 
cally simulate the configuration of Fig. [3^ by including the 
effect of a solid substrate below the sheet using a simple re- 
pulsive potential. Tearing is driven by a rectangular tool of 
width 15h that leans into the direction of motion at an angle of 
45°, just as in the experiment: this leads to stable concertina 
tearing as the tongue is confined between the tool and sub- 
strate. When the tool advances, we observe new folds form- 
ing in two steps: 1) the sheet ahead of the crack tip wrinkles, 
and 2) the crack tips advance by passing on either side of the 
wrinkle just formed, which now folds plastically and runs into 
folds formed in earlier cycles (Figs. (3}),c). The resisting force 
lf]~2l in a cycle of folding and tearing increases as wrinkles 
form, reaches a peak just before crack advances due to strong 
stretching, and then falls before the cycle starts anew; unlike 
in brittle tearing, fracture is quasistatic and tensile. This form 
of tearing is also insensitive on the form of tool used. 

To determine the wavelength of the folds as well as the 
crack paths, we must consider the stress field ahead of the 
two crack tips as a function of tool motion. Once a stack of 
folds has formed, the force due to tool motion is transmit- 
ted through these to the sheet ahead over a width w, which 
deforms out of the plane with amplitude A and wavelength 
A, measured in the direction of tool motion. The resulting 
stretching strain (A/w) 2 is primarily in the direction perpen- 
dicular to tool motion. The stretching energy in an area Xw 
scales as Us ~ Eh(A/w) 4 Xw while the bending energy 



in the same area Ub ~ Eh 3 (A/X 2 ) 2 Xw. Minimization of 
Us + Ub yields A w w^/h/A [16]. The wavelength, how- 
ever, is limited by A > 2 A which leads to an expression for 
wrinkle wavelength given by 

A = cxw^h 1 / 3 . (1) 



Further compression only folds the wrinkles plastically. To 
test the results of our theory and numerical simulations, we 
carried out tearing experiments using aluminum film of thick- 
ness h = 24 /im, yield stress a y « 0.003£? and toughness 
T « dyh 1 15 ]. Using image analysis, we extracted the wave- 
length A of the folds as a function of the width w of the tongue 
and compared them to the results of simulations (Fig. [3]i), and 
find an excellent match between simulation, theory and exper- 
iment, even including the weak dependence of A on the yield 
stress. Both experimental and simulation data are well de- 
scribed by equation ([I]) with the constant c\ = 2.1 ± 0.1 even 
for materials with a high yield stress when deformation ahead 
of the crack tips is primarily elastic. Interestingly, our elastic 
analysis for A leads to an expression that is identical to one 
based on a purely plastic analysis (HO, despite the different 
material assumptions, owing as is typical in these situations to 
the strong constraints imposed by the geometry of the system. 
For very thin sheets, when h/w —> our elastic analysis is 
always valid for since the maximum strain of bending during 
the wrinkle formation becomes vanishingly small compared 
to the yield strain. 

To understand the crack paths we assume, as is usual, that 
crack extension occurs when the energy available for crack 
growth is sufficient to overcome the resistance of the material 
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FIG. 4: (color online). Perforating (a) aluminum film or (b) simu- 
lated ductile sheet (cr y /E = 7 s /20 with 7 S = 0.4) displays radial 
cracking with four petals. Color indicates deformation energy den- 
sity. 

l6lfT71l. For our system this leads to the expressions 

^=2ThcosO (2) 
ox 

^ = -2rftsin0, (3) 
aw 

where equations (2) and (3) correspond to force balance paral- 
lel and perpendicular to the middle line, respectively, U is the 
elastic energy, and is the angle of crack paths with the mid- 
dle line. The factors of two follow from both cracks advanc- 
ing simultaneously. Dividing (3) by (2) we obtain tan# = 
dw/dx = — / f^- • Since the energy is distributed on an 
area of width A ahead of the crack tips, ^ « U/X, while 
|g = U/w so that ^ = tan (9 - A = C x(h/w)^ 3 9 where 
the second equality follows from eq. ([I]). Using w(0) = we 
can integrate the previous equation to obtain 

w = cohVW*. (4) 

This form with cq = 2.4 is in good agreement with both ex- 
periment and simulations (Fig. |3^). Once again, our result 
is analogous to the case of fully plastic tearing [3], where 
the above scaling relation arises by postulating that cracks 
advance perpendicular to tensile stress in a purely kinematic 
model. Again this is because geometry rules the crack path. 

We conclude with a discussion of tearing induced by the 
motion of a tool through the sheet in the out-of-plane direc- 
tion (Fig. |4^). This leads to petal cracks which has been stud- 
ied previously in ductile sheets using a combination of exper- 
iments and scaling concepts [QjOIIS). We complement these 
results using numerical simulations of tearing driven by a fric- 
tionless conical tool. The perforation force F = OUt/dR, 
where R is the radius of the intersection of the tool and orig- 
inal plane of the sheet and U t is the total energy of deforma- 
tion and fracture. As in concertina tearing, cracks are sup- 
plied with energy by stretching of the sheet between their tips. 
Assuming that the amplitude of the deformations is A, with 
A <C R implying a polygonal hole, the fracture energy is 
given by U c = nThR c , where R c = Rj cos 9 is the distance 
of the crack tips from the center of the perforating hole, with 
6 = ir/n the half angle between cracks. In the case of brittle 
fracture, the resisting force F w dU c /dR = nTh/ cos (tt/u). 
The number of cracks can be predicted by minimizing F 



which occurs for n = 4 and yields F = 5.66I7i. Again this 
result is to be contrasted with the ductile case where plastic 
bending of the petals contributes significantly to the force, and 
yields an optimum number of petals n = 4 lf]~8l suggesting a 
universal geometry consistent with our simulations for both 
brittle and ductile materials (Fig. j4j}). 

Our study highlights both the morphological complexity 
and the geometrical underpinnings of driven cracks in brit- 
tle and ductile thin sheets. Using numerical simulations, we 
have explored the forms that result and highlighted the qual- 
itative differences between the oscillatory tearing in brittle 
sheets mediated by a single crack, and the concertina folding 
and tearing mediated by two cracks in ductile materials. Sim- 
ple scaling laws allow us to explain our observations, which 
we also corroborate with experiments. Interestingly, our re- 
sults are similar to those obtained using purely plastic analyses 
and highlight the role of geometry again in these systems. A 
minimal phase diagram characterizes the morphological phase 
space, and suggests approaches for the controlled tearing of 
thin films for structure and function. 
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